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Abstract 

A general method for focusing of waves, based on phase and amplitude compensation, is applied 
to monochromatic, polychromatic and diffusive waves. Monochromatic waves may form spatially 
localized waves in free space, whereas polychromatic waves form non-decaying traveling evanescent 
modes confined to subwavelength regions in media where the frequency depends on the wave 
vector. We suggest an analogy between the phase compensation method and the transformation 
of frequencies between inertial, relativistic coordinate systems. 



I. INTRODUCTION 

Physicists and engineers often need to find methods for confining a field as much as pos- 
sible with a minimum of knowledge of the system at hand. In the case of monochromatic 
waves, it has been found that phase and amplitude wave engineering is of crucial impor- 
tance for subwavelength focusing of wavesp~H3]. In the case of pulsed waves, the situation is 
more complex, mainly due to the fact that the superposition of many spectral components 
in most cases tend to broaden the focal spot. Theoretically, one may hope to be able to 
focus strongly confined pulses as in Ref. [I], but the aperture field and medium required 
remains elusive. However, there are a few promising pathways, e.g., using rainbow-colored 
apertures which combine to a single white-colored focus j5] or diffractive optical elements[6]. 
Well developed theories and numerical algorithms are available for studying monochromatic 
and polychromatic waves in focal regions jH EH9], but they do not provide direct insight 
into the reverse engineering problem and may not provide a fast solution or be optimal 
with the minimal amount of information given. Of particular interest here is the emerging 
field of evanescent wave focusing, where up to now mainly monochromatic waves have been 
studied[3l fl0HT3] . However, also for waves containing several frequency components it is 
necessary to have straightforward methods which allow backtracking of the signal in a man- 
ner such that the initial field resulting in a strong focus is found. In the present work we 
detail a method for finding the intensity distribution near the focal region using phase and 
amplitude compensation. In order to achieve this goal, we pinpoint a certain distribution at 
focus and then try to calculate the required aperture field. Moreover, we show an analogy 
between our compensation method and the transformation of frequencies between inertial, 
relativistic coordinate systems. 

II. MONOCHROMATIC WAVES 

The method of phase and amplitude compensation has been described in a number of 
studies (see e.g. Ref . [H [2j [H] and references therein), but was more recently re- introduced 
to study focusing of evanescent waves [3J. The procedure is most easily demonstrated for 
scalar monochromatic waves satisfying the Helmholtz equation 
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where k = u/v is the wavenumber, u is the angular frequency and v is the velocity of light 
in the medium under consideration. Using the angular spectrum representation, the field 
can be written as 

E(x,z) = — / E k (k x y h * x+lk * z dk x , (2) 
where the angular spectrum is given by 

/oo 
E(x\0)e-^ x 'dx' , (3) 
-oo 



where 

' yjk 2 - h 2 x if k 2 > k 

k-, ~ 



iy/kl - k 2 if k 2 < k 2 £ 

Here k 2 > k 2 represent the homogenous plane waves, whereas k 2 < k 2 correspond to inho- 
mogeneous or evanescent plane waves that propagate in the x direction and decay in the z 
- direction. 

The trick is now to express the angular spectrum as 

E k (k x ) = f k (k x )e~ ik * z ° , (4) 

where we must again require E(k x ) to be bandlimited (i.e. nonzero only in a range k a < k x < 
kb) and converges sufficiently quickly such that the solution for the field does not possess 
any unphysical divergency problem. The conditions for this was discussed in Ref. [10], and 
will not be repeated here. The field at the focal point (z = z ) is now given by 

1 f°° 

E(x,z ) = — / f k (k x )e tk * x dk x , (5) 

i.e. a Fourier transform of the aperture function f(k x ) resulting in a finite resolution at the 
focal point. As an example, assume now that fk{k x ) = Eq for \k x \ < k a (k a k) and zero 
elsewhere. Then we have 

E(x,z ) = — / e lkx dk = E ^ '- . 6 

27T J_ kQ 7TX 

In order to find the field E(x, 0) we must evaluate 

i rko 

E(x,0) = — / f k (k x )e ik ^- ik ^dk x . (7) 

In Ref. [10] it was found that in the case of paraxial waves (k x k) the field E{x, 0) is 
approximately a quadratic phase function, which has been described in detail in standard 



textbooks on optics pQ. Thus, in the paraxial approximation the phase compensation method 
directly leads us to the well-known phase profile for a lens giving the most strongly focused 
profile[T]. Recently, the phase compensation method has been used to design evanescent 
waves in focal regions [3| [TO]. The simplest wave profile exhibiting a focal region is found by 
letting fk(k x ) = Eq if k a < k x < kb and zero elsewhere. Here k a 3> k and Eq is a constant. 
The field can then be expressed as 

E 



e ik b x-k b (z-z ) _ e ik a x~k a (z-z ) 



(8) 

IX — [z — z ) 

The intensity at the focus is therefore just ex sine 2 ( kb ^ ka x) whereas the aperture field 
oscillates rapidly within an envelope such that the intensity is cx l/(x 2 + z 2 ). Notice that the 
aperture field is much stronger than the field in the focal region, i.e. E(x, 0)/E(x, Zq) ~ e kbZ °. 
However, we also see that the intensity distribution is more confined at the focal region than 
anywhere else. In Fig. fiTa) the intensity of a monochromatic, evanescent wave is displayed 
with k a = 0.1 (arbitrary units), kb = 1 (arbitrary units) and k^ZQ = 5. The solid line shows 
the field in focus, whereas the dashed line shows the fields at z = 0. This illustrates our point, 
namely that we may obtain a more spatially confined intensity distribution (although the 
magnitude is considerably smaller than that at the aperture). If we assume that kb 3> k a , a 
possible criterium for focusing could be obtained by requiring that the aperture intensity (at 
z = Zq) envelope at x = z , corresponding to the position with half the maximum intensity, 
should be wider than the focused intensity distribution of half-width ~ Tr/kb- Thus, by 
requiring kbZo > 7T, we ensure that the field is focused. At the same time, the field should 
not diverge at the aperture, and it is clear that fc&zo cannot be too large. 



III. DIFFUSIVE WAVES 



Interestingly, the compensation method described above is not limited to designing waves 
in focal regions described by the Helmholtz equation. It can also be used to construct so- 
lutions to other differential equations, such as the diffusion equation, and we will here look 
at one example. Imagine that we generate a spatial magnetic field distribution by position- 
ing a system of co-aligned current-carrying wires in a homogeneous medium of frequency- 
independent conductivity a and permittivity e. Each of the wires carry a low frequency 
current where the frequencies of any temporal waveform fulfills u> <C a/e. In absence of 
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sources the resulting electric (or magnetic) field can be described by the diffusion equation 
on the formf 



fl^,.)^. (9) 

where D = l/(ere) is the diffusion coefficient. We now want to describe the field in the 
vicinity of the wires, but not including any sources such that Eq. [9] is valid. The field can 
then be written on the following general form 

E(x,t) = — E k (k)e lkx - Dk2t dk , (10) 

2n ./_™ 



where Ek(k) is the fourier amplitude. Usually, a solution of eq. 10 will represent a diffusion 
process where the field delocalizes with increasing time [18]. However, we are here interested 
in finding a field E(x, 0) with an optimally localized field E(x, to) {to > 0). To see how one 
may design such an initial field, let us assume that E k {k) = E exp{Dk 2 t ) is a constant for 
| A; | < fco and zero for \k\ > ko- Thus, the field is given by 



E(x,t) = ^ / e tkx ~ Dk ^dk , (11) 

27r -/-fco 



The real value of eq. 11 at t — is given by 



2vr 

whereas at t — to it is 



E( x ,0) = ^ I cos{kx)e Dk2to dk , (12) 

fc 



E(x,t ) = ^ [ k °e ikx dk = Eo Sin{koX) , (13) 



2tt J_ kQ TlX 

By selecting Dto such that Dk^to ^> 1, we see that the leading contribution to E(x,0) ~ 
cos(kox)e Dk o to , which is an oscillatory function with period Ax = 2ir/k . Thus, the field at 
t = is strongly delocalized, but with sharp local fluctuations, as seen in Fig. [I]b). On the 
other hand, E(x,to) is strongly localized with a central lobe lixjko-, but is a factor e Dk o to 
smaller than the field at t — 0. A sharp localization is therefore obtained at the expense 
of an exponential reduction in field amplitude. This is rather similar to the behavior we 
observe for monochromatic evanescent waves, although the details are different. In addition 
to describing electromagnetic fields, the diffusion equation may also describe concentration 
or temperature (but now with a different diffusion coefficient). In those cases it should be 
mentioned that we do not consider absolute temperature or concentrations (as required by 
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our boundary conditions), but rather modulations about an equilibrium. Thus, the negative 
values seen in Fig. [l]b) do not mean absolute negative temperature and concentration, they 
are just signatures of the variations about the mean temperature or concentration seen when 
designing an initial profile required to bring the diffusive waves to a focus. 



IV. POLYCHROMATIC WAVES 



An interesting question is whether the method above can be applied directly to polychro- 
matic waves (e.g. pulsed waves) satisfying the wave equation 

^9 x 1 d 2 E(x, z,t) , , 

V 2 E(x,z,t) = - (14) 

In the angular spectrum representation, the field can be expressed as 

poo POO 

E(x,y,z,t)= / / E k (k x ,u)e ik ' x+ik ''- iu ^dk x du , (15) 

JO J-oo 

where k = u{k x )/c and 

' y/k 2 - kl if k 2 > k 2 
i^k 2 x - P if k 2 < k 2 x 

In order to bring the polychromatic waves to focus we require cross-compensation of the 
phase, which amounts to setting k z z\ — uj{k x )ti = 0, i.e. the spatial part of the phase is 
compensated by the temporal part at (x,zi,ti). Such a requirement can only be fulfilled if 
u depends on the spatial frequencies (i.e. direction of each plane wave), where the angular 
frequency is given by u(k x ) = i^{yi)\k x \vi. Here 7(1*1) = — [vi/c) 2 and v\ = z\jt\. 

The fact that the frequency of the spatial wave vector depends on the direction was utilized 
in Ref. [S] to combine a rainbow spectrum to a single, focused spot. Conceptually the 
idea presented here is somewhat similar, although we use an entirely different approach to 
achieve the goal. Note that we must distinguish between the two cases v 1 > c and V\ < c. 

The case v± > c results in superluminal, localized waveforms which have been discussed 
extensively in the literature [T6l [17]. Here oo(k x ) is real and given by u(k x ) = ak x (a = 
v i I \j (vi/c) 2 — 1) . The scalar field is then given by 

/•OO ^OO / \ 

E(x,z,t)= / E k {k x ,u)e ikxX+ilkxla V-^)dk x du . (16) 

JO J -00 
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A similar representation was found in ref. [T7] describing one-dimensional nondiffracting 
pulses, but here we have arrived at Eq. [16] from an entirely different perspective based on 
requirement of specific field profile resulting from phase compensation. 

Of greater interest in the current study is the case v\ < c, since that represents an 
evanescent wave solution similar to that seen in the previous section. Now a is purely 
imaginary, such that 

/*oo /*oo / \ 

E(x,z,t) = / E k (k x ,u)e ikxX ~ lkxh V^)dk x du. (17) 

J0 J -co 

As an example, we may approximate the angular spectrum as E(k x ,u) = Eq for k a < k x < 
kb , < u < u)2 and zero elsewhere. Here we also assume that k a 3> k, such that only 
evanescent waves are excited with k z ~ i\k x \ and u(k x ) ~ i\k x \vi. This evanescent-wave 



approximation therefore requires v i <C c, and the field will be 

' ^ik b x-k h {z-v\t) ( J,k a x-k a (z-v 1 t) 

E(x, z,t) oc 



ix — (z — fit) 



(18) 



Note that eq. 18 is identical to eq|8]if we set Zq = V\t, and is therefore just a translation of 
the evanescent wave along the optical axis with the intensity distribution similar to that of 
fig. [!} Since we require k b k, which is fulfilled if we, as an example, set kb = 10k = 20n/X, 
the half-width of the intensity distribution at the traveling focal line is narrower than 
~ n/kb = A/20. The intensity at z = grows as I(x, 0,t > 0) oc exp(kbVit) / (x 2 + (fit) 2 ). 
However, when t ^> 1/kbV ~ t\ our approximate theory for evanescent waves does not longer 
hold in the vicinity of z = (this follows from the considerations above; see also Refs. [TUirH] 
for a discussion about diverging solutions). The fact that the evanescent wave does not 
change at the focal line z\ = v\ti is surprising, given the condition E(x,0,0) = E(x,Zi,ti) 
above. However, we also note that E(x, 0, t±) 3> E(x, z±, ti) such that the increasing energy 
of the field at the aperture is used to keep the field at the focal line z\ = viti unchanged in 
magnitude as it propagates outwards. The technical implementation of focusing of polychro- 
matic waves is in effect similar to that of monochromatic waves, but with two new important 
features: a) The focusing must take place in a spatially dispersive media, and b) The field 
at the aperture must grow with time (within the approximation given above). 

An interesting analogy occurs if one compares the problem of designing an aperture field 
that generate optimal focus with that of transformation of frequencies between inertial, 
relativistic coordinate systems [IS]. To see this, consider a plane wave of the form exp(k z z\ — 



u(k x )ti + <p(k x )). Any wavepacket is a weighted sum over such plane waves in spatial and 
frequency coordinates. In order to bring the wave packet to focus we require that every 
spectral component is at focus a certain distance Z\ from the aperture at a given time 
t\. This can be done by setting k z z± — u(k x )ti + <p{k x ) = 0, such that the phase is exactly 
compensated at (z\, ti). It should be noted, as seen above, that such a requirement may give 
rise to strong localization at other positions as well. In any case, such a cross-compensation 
requirement can only be fulfilled if u depends on the spatial frequencies (i.e. direction of 
each plane wave), where the frequency is given by 



uj{k x 



<f){k x )/t 



i ± 4>(k x )/h 



(-) 2 "1 (-) 2 "l\ 



1 + 



Vl 



1 



k x zi 
4>{k x ) 



1 



(19) 



Now consider the problem of transforming plane waves between inertial, relativistic coordi- 
nate systems. That is, consider an observer moving at a speed v relative to a fixed frame. 
According to Einstein's special theory of relativity, the moving observer will detect a fre- 
quency {J given by tu' = 7(f) (u — k z v), where 7(f) = 1/ 



T — (f/c) 2 [T3]. From the above 
it may be inferred that this expression corresponds to k z Z\ — u(k x )ti + i<p{k x ) = if we 
make the associations 4>{k x )/ti — > u}'/j(vi) and v\ — > v, i.e. the phase factor 4>{k x ) must be 
associated with the frequency measured in the moving frame. In the special case considered 
in this study (j)(k x ) = 0, which corresponds to w' = 0. For evanescent waves we then imme- 
diately find the spatial dispersion relationship given above. For propagating waves we may 
set k z = (u/c)cos6, where 9 is the real angle at which a specific plane wave makes with the 
direction of motion. Then we must have cos9 = c/v and therefore v > c, in agreement with 
the observations made above. 



V. CONCLUSION 



In conclusion, we have suggested a compensation method for designing aperture fields 
giving rise to strongly confined waves. The idea is to first try to compensate the phase 
or amplitude such that only transverse spatial frequencies are left in the angular spectrum 
representation at a given region in space, thus allowing a strongly focused wave to form here. 
Next we calculate the aperture field required to obtain such a phase or amplitude compensa- 
tion. The method has been employed to study diffusive, monochromatic and polychromatic 
waves, and shown to give new insight into the problems at hand. The method here can 
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probably also be applied to other wave systems where phase or amplitude compensation is 
beneficial. 
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FIG. 1: In (a) the intensity of a monochromatic, evanescent waves is displayed with k a = 0.1 
(arbitrary units), k^ = 1 (arbitrary units) and k\,ZQ = 5. In (b) the diffusive waves (a) with ko = 1 
(arbitrary units), D = 1 (arbitrary units) and to = 10 are displayed. The solid lines show the 
field in focus, whereas the dashed lines show the fields at t = for diffusive waves and z = for 
monochromatic evanescent waves. 
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